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1. INTRODUCTION 

Nowadays, integrated circuits electromagnetic susceptibility threshold gets lower due to 
the integration densities increasing and higher levels of metallization. This increased vulnerability results 
from their size shrinking, supply voltage and the regular increase in their operation frequencies [1, 2]. 
Gigahertz frequencies are particularly harmful for these electronic systems because their wavelengths are 
likely to generate resonance phenomena on the integrated circuits tracks, thus increasing the system 
disturbance risks [3-5]. 

As microelectronics technology progresses, manufacturers are continually striving to build higher 
density electronic systems integration. This evolution is characterized by a regular device shrinking and 
a multilevel integration [6]. This led to acomplex interconnection circuitry consisting of parallel planes of 
equidistant transmission lines. Today this becomes a complex problem faced in the improvement of 
the integrated circuits performances [7]. When designing and implementing an electronic system, 
the possible electromagnetic interactions should be taken into account and methods should be provided to 
reduce their influence sufficiently, so as to ensure safe operation in most cases [8, 9]. 
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Ensuring electromagnetic compatibility at the integrated circuit level implies effective reduction of 
noise sources and disturbances origins. The modeling of these facts becomes mandatory to ensure 
a predictive approach to the risks of EM-induced disruption [10, 11]. This requires specific tools, models and 
EMC knowledge [12, 13]. The designer of integrated circuits should consider the capacitances and 
inductances that were not considered in earlier times, because of the interconnection’s lengths and the higher 
operating frequencies [14, 15]. The study of electromagnetic radiation in an integrated circuit requires 
knowledge of the electromagnetic fields in this device and therefore the flux values [16]. The calculation of 
the magnetic field requires a resolution of Maxwell equations. This is achieved by several computational 
methods, such as the Method of Finite Element, finite differences and the method of moments [17, 18]. 

In order to further improve the integrated circuit several calculation models have been proposed to 
estimate the behavior of the interconnection’s lines, as it was presented in [19. 20]. And may researchers 
have calculated and simulated the electromagnetic problem [21-23]. In this paper we have used the finite 
element method to determine the radiation and the coupling between its interconnection tracks in order to 
integrate all the line fields surrounding an interconnect line. We have organized the paper as follows: in 
the seconde part we have given the definition of the finite element method, in the third part a mathematical 
model has been developped to minimize electromagnetic disturbance on circuits interconnects, and in 
the fifth part the circuits interconnect has been simulated with defferents parameters in order to achieve 
an interconnection structure with less energy losses. 


2. THE FINITE ELEMENT METHOD 

The finite element approach is based on solving partial differential equations knowing the boundary 
conditions. This method was initially used for solving problems in the field of fracture mechanics and 
structural design (by mechanics). The method was first used for computing the electromagnetic field in 
the 1970s by P.P Silvester and M.V.K Chari. In most cases, it 1s integrated with C.A.O softwares and showed 
great advantage for the designers of physical systems [24]. 

The basic approach of the finite element method is to subdivide the field of study into finite 
numbers of subdomains called elements as show in Figure 1. The approximation of the unknown values is 
done for each element of the interpolation functions. This function is defined based on the geometry of 
the element that is chosen beforehand and fits with the nodes of this element with respect to the unknown 
values. This is called nodal interpolation [25]. 
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Figure 1. Mesh by finite element 


3. THE MATHEMATICAL MODEL 

In the following, we will be developing a mathematical model that governs electromagnetic 
compatibility problems of the magnetization conservation. Using the Maxwell equation, we can deduce 
the relation that links the magnetic induction to the magnetic vector potential as: 


V.B=0>454/B=V*A (1) 


introducing (1) into Maxwell Faraday equation we will have: 


=>. 


(245) ag (2) 


since the rotational of a gradient is null, we can deduce the existence of an electric scalar potential V, and 
then we write it as: 
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IV/E + — = —gradV (3) 
Ot 
we have: 
F (4) 
ðA 





E = -—-— V 
AL grad 


replacing (4) in the expression of Ohm's law we get: 


> 


P eee ðA 

J =ogradV — oF, (5) 
we have 

J. = ogradV (6) 


J: current density of the excitation source 








> > ð 
=J,—-o— (7) 
J Je o Ot 
by introducing the (1) in the relation of magnetic medium we get: 
és VA 
H = (8) 
H 
replacing (8) in the Maxwell-Ampere equation and applying to it the rotational we will have: 
Be N'A > ð aD 
y^ =J,-o—+—— i 
( 7 ) E (9) 


taking into account the relation of a dielectric medium and replacing (4) in (9) we obtain the following set 
of equations: 
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it is known that: 


V^V^A = -V2A + grad (V. A) 





(13) 
E E. 32A E 
-IPAR -oikene grad (VA + we) (14) 
to completely define the magnetic potential vector, one imposes the condition of Lorentz: 
O OV 
V.A + ue — = 0 (15) 
rE at 


Then (14) is simplified as given below: 
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Logs 0A FA 
u Jt apn 


(16) 


we consider that the propagation is done in a Cartesian coordinates system along the z axis, and then 


the propagation equation of the magnetic vector potential will be: 


u 
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if we consider the sinusoidal harmonic regime this leads to: 


1A, Az) ve 
u Ox2 dy? Gow EW ) Z = J, 


(17) 


(18) 


Az, Jz: are respectively the component of the magnetic vector potential and the current density source along 
the z axis. In the (18) characterizes the propagation of the magnetic vector potential. This latter will be 
considered in the following part of this work to model the electromagnetic interaction phenomena of 


the interconnections in an integrated circuit. 


For solving the potential magnetic vector propagation equation, considering the finite element 
formulation, we have introduced the Galerkin method [14] defined in a Cartesian coordinate system. 
The propagation is considered according to the oz axis, by considering the phenomenon is governed by 


the equation: 








1 ae 32A, 


wl ax2 ` ay? 
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I) (- rfe Gye + Joo- ev Jz) = Pi dx dy = 0 
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|+ (jow — ew*)A, — J, = 0 





using the boundary Dirichlet conditions, this leads to: 
| oA: ar =0 
BZ on Hg 


-7 || PoAzdxay + (jow — £w?) || Actiaxay = || Jedidxay 


N 


A= X bA 


j=l 


3 -4 f| 96.96 axay| A; + Gow — a) ff #6 ,axay| A; = || J-bidxay 


for all mesh nodes, the following matrix system can be written: 


[M][A] + Gow — ew*)[L][A] = [K] 


1 — — 
HJJo 
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Lij = If pjdxdy = || diteaxay (29) 


[A] and [K] are respectively the unknown vector and the excitation source vector 
[M] and [L] are two matrices that depend on the mesh structure 


4. PARAMETERS IDENTIFICATION OF THE EQUIVALENT CIRCUIT DIAGRAM 

It is important for an integrated circuit designer to predict the behavior of his product in an 
electromagnetic interference environment. Interconnections are often victims of this electromagnetic 
pollution. This part consists of a numerical modeling of the electromagnetic interferences. We will analyze 
the influence of the frequency parameters of the supply signal and spacing between track and the permittivity 
of dielectric on the evolution of the parameters of the equivalent diagram. The interconnection lines can be 
modeled using their distributed constants, the equivalent diagrams are shown in Figure 2, it represents a track 
portion in the form of an electrical quadrupole that links the characteristics (1, v) of the position (z + dz) to 
that of position (z) at time (t). 


(+ d > R 


—— p Ved ee 
v(z,t) v(z+dz, t) : 


v(z,t):represents the voltage at position "z" 
v(z+dz,t):represents the voltage at position "z+dz" 


Figure 2. Equivalent circuit of a strip portion 


The parameters of the equivalent diagram are: the equivalent inductance "L" per unit length, it 
characterizes the magnetic energy density stored in the medium, the equivalent capacitance "C" per unit 
length, it characterizes the density of dielectric energy stored in the substrate, the resistance series per unit 
length "R", it characterizes the losses by Joule effect, the parallel conductance per unit length "G", it 
characterizes the losses in the insulation. Joules losses were used to evaluate the linear resistance given by 
the following relation: 


P 

R= > [2/m] (30) 
T 

LE mas [H/m] (32) 


to evaluate the inductance, the following magnetic energy should be used: 


1 
Wnag = [|| ZrorlHlèav (33) 
2* Wai 
=— 5 [F/m] (34) 


the dielectric energy is used to evaluate the linear capacitance: 


Waietc = ||| eoerlel2av (35) 
v 


to evaluate the conductance of the substrate, we will use a formula that includes the capacity: 
G=C*tw*tand  [S/m] (36) 


tano: dissipation factor (dielectric loss angle) 
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To study the parameters of equivalent diagrams as a function of the supply frequency, the distance 
between the track and the permittivity, we have elaborated the flowshrt shown in Figure 3 which has been 
used to design our computation programs. At the start of our program as shown in Figure 3 we introduce 
geometry and physical characteristics from model to modeled. After having meshed the field of study, we 
will solve the partial differential equation at each mesh point and identify the parameters (R, L, C, and G) of 
the interconnection line for a given frequency and distance between tracks. Each time the frequency and 
the distance between tracks are incremented. 


Introduction of geometry 


Introduction of physical 
characteristics 


Incrementation 
Of the distance 


Between strips 
Resolution of the partial derivative 
Equation in terms of the magnetic vector 
Incrementation 
offrequency 


Identification of the parameter of the 
equivalent schemas of the microstrip 
line (R, L, C, G) 








NO 
End of change 
the distance 
beteen strips 


Yes 


Figure 3. Flowchart used to study the parameters of the equivalent circuit as a function 
of supply frequency, the distance between the track and the permittivity 


5. RESULTS AND DISCUSSION 
5.1. Ground plane effect 

In this part we will study the effect of ground plane on the impedance of the interconnects. 
We consider the structures in Figure 4, this microstrip coupled interconnect has the following geometrical 
parameters w = 2 um, S = 5 nm, h = 1 um, H = 3 um, e = 3.9 In this case we note a significant increase 
in the capacity towards the mass which is due to the lines of fields which refocus towards the mass as show 
in Figure 5. 


5.2. Influence of the frequency 

In our analysis we based our system on a two parallel transmission lines. This corresponds to a 
2 lines of an n line internal integrated circuit bus. The three-line configuration is used to discuss 
the interdependance of the two-line Transmission lines. This microstrip coupled interconnects have 
the following geometrical parameters: w = 2 um, S = 5 nm, h = 1 um, H =3 um, e£ = 3.9. The three microstrip 
lines as shown in Figure 6 are equidistance and mounted on a oxide grown on silicon substrate. We simulated 
and modeled the interconnect capacitance in different situation. The Figures 7 show inductance and 
resistance versus frequency.The increase in resistance is due to the skin effect which causes the distribution 
of the current density in the section of the track to become non-uniform with the increase of the frequency 
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and the current that travels the track will tend to circulate on this periphery in a section defined by a thickness 
called "skin thickness". Therefore, the active section of the track decreases and therefore the resistance 
increases. 

We also note that the inductance of the track decreases as the frequency increases as shown in 
Figure 8. According to the relation of the inductance it is understood that, as regards the internal volume of 
the conductor, the contribution to the integral of the magnetic energy is even smaller than the frequency is 
large. It is found that the capacity increases with the increase of the frequency, this is explained as much as 
the frequency is more important the dielectric energy will tend increased consequently the capacity increases. 

The increase in conductance as show in Figure 9 is due to the increase in capacity, According to 
the equation of conductance. The obtained results well motivate the effect of the electromagnetic interactions 
on the electrical parameters of the interconnects, notably the impedance which has known an increase 
compared to that obtained in the absence of the electromagnetic interactions as show in Figure 10. 


x 10° 








Capacitance parasitic[pF/m] 








(a) (b) 
frequency [Hz] 
Figure 4. Microstrip coupled interconnect: Figure 5. Parasitic capacitance 
(a) without ground plane, (b) with ground plane versus frequency 





ENNANANAAAANANAANAAAAAARAAARRAAARRARA RRA RRA RRS RRA SARRRNE EELS 


Figure 6. Symmetric microstrip coupled interconnects 
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Figure 7. Inductance and resistance Figure 8. Capacitance 
versus frequency versus frequency 


Investigation on EM radiations from interconnects in integrated circuits (Lounas Belhimer) 


308 O ISSN: 1693-6930 


250 1 


ni '|-0—] strips 
—Ə— 2 strips 
—O— 3strips 


200 


an 
C 
Impedance 
(m) 
oo 


= 
Oo 
Oo 


conductance[S/m] 


50 





104 108 108 1010 10° 10? 104 10 108 1010 
frequency [Hz] frequency [Hz] 


Figure 9. Conductance versus frequency Figure 10. Summary and data comparison 


5.3. Influence of the dielectric permittivity 

We varied the dielectric permittivity and the frequency is fixed we calculate the new values of 
the parasitic capacitances. The permittivity is directly related to the dielectric energy so by increasing 
the permittivity the dielectric energy also increases consequently the parasitic capacitance increases as show 
in Figure 11 and Figure 12. 


-15 „40-15 
6 x10 . 7, 10 


p [Si] © 


Ww 


Capacitance parasitic[pF/m] 
Capacitance parasitic[pF/m] 





Permittivity[F /m] Permittivity[F/m] 
Figure 11. Parasitic capacitance Figure 12. Parasitic capacitance 
versus Permittivity versus permittivity for different S value 


5.4. Influence of the distance inter tracks 

For different values of "S" (the distance between two interconnection tracks), it sets the working 
frequency is | gigahetrz we calculate the value of the capacity and inductance per unit length. As show in 
the Figure 13 and Figure 14, the capacity decreases with the distance between tracks, this evolution can be 
explained by the fact that we associate the increase in the distance between tracks, in order to locate an 
optimal layer thickness corresponding to the thickness for which the spurious capacitance track ratio is 
the lowest. 


5.5. Modeling of crosstalk 

Any unwanted voltage or current and created by coupling between two tracks of an interconnection 
network is considered to be crosstalk. In what follows, we propose to model the tension induced by a track on 
another track that is adjacent to it. For that we consider the precinct structure. According to the law of 
faraday, the coupling in magnetic field induces a voltage Vin directly related to the variations of the magnetic 
flux ® which crosses the conductor of section "s" and it is written as follows: 


dp 
Vin = Ge (37) 
$ = || Bas (38) 
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B=rotA (39) 


In this part always considers the Figure 6, but this time with the following dimensions: w = 20 um, 
S = 20 um, h = 5 um, H = 40 um, e = 3.9. Considering a supply voltage of 2 V, with a track which has a 
length of 0.5 m, we get the following simulation results in Figure 15 and Figure 16. The Figure 15 represents 
the results of simulation of the crosstalk as a function of the frequency of the signal compared with that 
obtained in [4]. We see that there is some agreement between the two results, and a relationship of 
proportionality links the crosstalk and the frequency. We see that from the frequency 10 Hz the voltage drop 
begins to be felt. So, the signal injected at the beginning of the track undergoes attenuation due to the effects 
of high frequencies and electromagnetic interactions. Figure 16 represents the simulation of the voltage drop 
in the interconnection line for frequencies given as a function of the distance between tracks, it’s noted 
that the voltage drop decreases with the increase in the distance between track and it increases with 
increased frequency 
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Figure 15. Crosstalk versus frequency Figure 16. Pressure drop versus interline distance 


6. CONCLUSION 

The main conclusions are drawn for the study of the influence of dielectric width and permittivity. 
The transmission line length and geometry was also studied and the impedance per unit length is treated in 
order to estimate the losses and define the limit amplitude of the digital signal in use. We demonstrated that 
the capacitance and resistance increase with increasing frequency but the inductance decrease and by 
increasing the permittivity or the distance between tracks the parasitic capacitance decrease. Obtained 
literature results and with our simulations data, demonstrated that the crosstalk as a function of the frequency 
of the signal. Through these results, we can estimate the best materials for the best circuits interconnect 
geometry to ensure the best possible performance. The results are found to comply with literature and 
encouraging for further investigations on dielectric separation. 
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